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Abstract. We determine the range of neutrino masses and cosmic radiation content 
allowed by the most recent CMB and large-scale structure data. In contrast to other 
recent works, we vary these parameters simultaneously and provide likelihood contours 
in the two-dimensional parameter space of Ncs, the usual effective number of neutrino 
species measuring the radiation density, and '^m^. The allowed range of '^rui, and 
Ncs has shrunk significantly compared to previous studies. The previous degeneracy 
between these parameters has disappeared, largely thanks to the baryon acoustic 
oscillation data. The likelihood contours differ significantly if ^ resides in a single 
species instead of the standard case of being equally distributed among all flavors. For 
^ = we find 2.7 < N^s < 4.6 at 95% CL while J^^-u < 0.62 eV at 95% CL for 
the standard radiation content. 
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1. Introduction 

The recent release of the 3-year WMAP data has stimulated several renewed analyses of 
cosmological neutrino mass limits. Neutrinos are known to have mass from oscillation 
experiments so that the unknown overall mass scale is unavoidable as a fit parameter of 
the standard cosmological model. The resulting mass hmits range from < 2.0 eV 

(95% CL) using the WMAP-3 data alone [1] to T.m„ < 0.17-0.4 eV (95% CL) when 
data from the Lyman-a forest is included [2,3]. 

Translating cosmological limits on the hot dark matter fraction into neutrino mass 
limits depends on the cosmic neutrino density that is fixed by standard physics and thus 
not an ordinary cosmic fit parameter. On the other hand, the most direct evidence for 
the presence of the cosmic neutrino sea derives from big-bang nucleosynthesis and from 
cosmological parameter fitting so that it is a natural consistency test to study if the 
cosmic radiation density implied by the cosmological precision parameters reproduces 
the standard radiation content. 

However, since neutrinos are known to have mass, one can not simply assume that 
J2'^u = when extracting an allowed range for iVcff. the effective number of neutrino 
species that is the usual measure of the radiation content. Nevertheless this has been the 
standard procedure in most of the recent parameter analyses based on WMAP-3 [2,4]. 
The caveat is especially relevant because wc found a degeneracy between Yl f^v = and 
A^cff) based on the cosmological data available in 2003 [5-8]. One result of our present 
J2 m,y-A^efr-analysis will be that this degeneracy is no longer present in the much smaller 
allowed range based on the 2006 data (Fig. 2). We will thus conclude that the current 
cosmological data provide essentially independent limits on X^m,^ and A^eflf- 

Another important issue is what one actually means with A^eff, as there are several 
different plausible cases that should be considered. One possibility is that the cosmic 
number density of the standard neutrinos is different from what is usually assumed, 
i.e. Y.f^v is equally distributed among all species that comprise A'cfr (our Case 1, see 
Table 1). Our second case is that three standard massive neutrinos have equal masses, 
i.e. the number of equally massive species is = 3 so that N^^ — N^ signifies additional 
radiation in some completely new form unrelated to ordinary neutrinos. Finally, we 
consider A^m = 1 that could represent a situation where the standard neutrinos are 
nearly massless, i.e. they have hierarchical masses with a largest mass eigenvalue given 
by the atmospheric scale of about 50 meV, while there is an additional massive species. 



Case 


Model 


Em^ (95% CL) 


A^eff (95% CL) 


1 


A^m = A^eff 


< 0.62 eV 


2.7 < A^eff < 4.6 


2 


iV^ = 3 


< 0.57 cV 


3.0 < A^cff < 4.6 


3 




< 0.41 eV 


2.7 < N^s < 4.6 



Table 1. Cases of neutrino mass distribution among the A^eff species. For the allowed 
E 'mv range we have marginalized over Ngs and vice versa. 
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perhaps a sterile neutrino. This case is largely motivated to demonstrate that the current 
cosmological data are sensitive to the Yl f^u distribution among the flavors. 

We begin in Sec. 2 with a brief description of the cosmological data used in our 
study. In Sec. 3 we derive the range of and A^eff allowed by these data and conclude 
in Sec. 4 with a summary of our findings. 

2. Cosmological data and likelihood analysis 

In order to study bounds on Y^rUi, and iVeff we use the same data as in Ref. [3]. We 
use distant type la supernovae measured by the SuperNova Legacy Survey (SNLS) [9] 
and large-scale structure data from the 2dF [10] and SDSS [11,12] surveys. From the 
SDSS we also include the recent measurement of the baryon acoustic oscillation feature 
in the 2-point correlation function [13]. Finally, we include the precision measurements 
of the cosmic microwave background anisotropy from the WMAP experiment [4, 14, 15], 
as well as the smaller-scale measurement by the BOOMERANG experiment [16-18]. 

We do not include data from the Lyman-a forest in our analysis. These data were 
used previously and very strong separate bounds on and A^eff were obtained [2]. 

However, the strength of these bounds is mainly related to the fact that the Lyman-ct 
analysis used in Ref. [2] leads to a much higher normalisation of the small-scale power 
spectrum than the WMAP data. Other analyses of the same SDSS Lyman-a data find 
a lower normalisation, in better agreement with the WMAP result [19-21]. In this 
case the Lyman-a data add little to the strength of the neutrino mass bound [3]. The 
discrepancy between different analyses of the same data probably points to unresolved 
systematic issues so that we prefer to exclude the Lyman-o; data entirely. 

We then perform a likelihood analysis based on a flat, dark-energy dominated 
model characterised by the matter density the baryon density 0^, the dark energy 
equation of state the Hubble parameter Hq, the spectral index of the primordial 



parameter 


prior 




Vl — Q^jn + ^DE + ^1/ 


1 


Fixed 




- 1 


Top hat 


h 


0.5 - 1.0 


Top hat 




0.014 - 0.040 


Top hat 


wde 


-2.5 - -0.5 


Top hat 


ris 


0.6 - 1.4 


Top hat 


as 


-0.5 - 0.5 


Top Hat 


T 


- 1 


Top hat 


Q 




Free 


b 




Free 



Table 2. Priors on the parameters used in our likelihood analysis. 
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power spectrum n^, the running of the primordial spectral index a.,, and the optical 
depth to reionization r. Finally, the normalization of the CMB data Q and the bias 
parameter h are used as free parameters. The dark-energy density is given by the flatness 
condition f^DE = 1 — — Including the neutrino mass Y^f^vi parameterised in 
terms of the contribution to the present energy density VL^h? = Y^rUy eV, and the 
effective number of neutrino species N^^, our benchmark model has 11 free parameters. 

Our priors on these parameters are shown in Table 2. The treatment of data is 
exactly the same as in Ref. [3]. When calculating constraints, the likelihood function 
is found by minimizing over all parameters not appearing in the fit, i.e. over all 
parameters other than and A^^efr . 




Figure 1. The 68%, 95%, and 99% confidence level contours for our three cases. 
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Figure 2. The 68% and 95% confidence level contours for Case 1 overlayed on the 
same contours using the data available in 2003 [6]. 



3. Bounds on neutrino properties 

Following these procedures we find the 68%, 95%, and 99% likelihood contours for 
and A^eflf shown in Fig. 1 for the three cases discussed in the introduction and shown 
in Table 1. The top panel of Fig. 1 corresponds to a nonstandard number density of 
the standard neutrinos, assuming a standard velocity dispersion as in all other cases 
as well. In Fig. 2 we overlay these contours with the analogous ones that we found on 
the basis of the data available in 2003 [6]. The allowed range of both parameters has 
shrunk dramatically as expected. Moreover, the pronounced degeneracy between X] 
and iVeff that was present at that time has now completely disappeared, largely thanks 
to the baryon acoustic oscillation (BAO) measurements. We conclude that at the level 
of precision that has now been reached, the cosmological data constrain and iVeff 

almost independently of each other. 

Perhaps the physically best motivated case is No. 2 where we have the ordinary 
neutrinos with mass (A^m = 3) and additional radiation in some new form. In this case 
we have a hard lower limit iVeff > = 3. Otherwise the contours of the middle panel 
of Fig. 1 are very similar to Case 1 (top panel) . 

The largest modification appears in Case 3 where we assume that all hot-dark 
matter mass resides in a single neutrino species. The mass limits are significantly more 
restrictive in this case. The reason is that for a single massive species the total neutrino 
energy density is larger in the semi-relativistic regime than if the mass is shared between 
all flavours [22]. Since the mass bound is such that neutrinos become nonrelativistic 
very close to the epoch of matter-radiation equality, this equality occurs later in the 
model with only one massive neutrino. As a consequence, small-scale structure is more 
suppressed, but the effect can be offset by a slight increase in the matter density. We 
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indeed observe that the best-fit value of 0^ is higher for = 1. However, both 
the SN la and BAO data prefer a low vahie of fl^n ^-nd consequently the model with 
iVm = 1 becomes a poor fit to this data at J2mi, around 0.4-0.5 eV. Table 3 shows 
exactly this effect. Here, Ax^ for Y^rrii, = 0.45 eV has been broken down into individual 
contributions from the different data sets. As expected, the main effect comes from 
SN la and BAO data. 



Data set 


Nm = N,s 


Nm = l 


CMB 


-1.2 


-0.9 


LSS 


0.5 


0.5 


SN la 


1.0 


2.3 


BAO 


1.3 


2.9 



Table 3. Ax^ for O^/i^ = 0.005 compared with the best fit model, broken down into 
individual contributions. 

We finally note that for the case with one sterile massive state and three active, 
almost massless neutrinos (the LSND 3+1 case) the mass bound is 0.45 eV at 95% CL 
(0.93 eV at 99.99% C.L.), somewhat lower than the 0.62 eV bound in the standard case. 
That the bound on the 3+1 model is stronger than for the standard case is contrary 
to what was found in previous studies (see [5-7]. The reason is the low value of 
preferred by the BAO and SNI-a data. 

4. Discussion 

We have derived likelihood contours in the two-dimensional parameter space spanned 
by a-nd A^eff, based on the latest cosmological precision data, however excluding 

Lyman-CK. We consider two physically motivated cases for A^eflf where either the effective 
number of massive neutrinos differs from the standard scenario, or where there is a new 
form of radiation besides N^ = 3 standard massive neutrinos. The results for these 
cases differ very little, except that in Case 2 there is hard lower limit A^eff ^ A^m = 3. 

For the sake of principle we have also considered a third case where all the neutrino 
mass resides in a single species. Here, the mass limit is more restrictive, refiecting that 
near the limiting mass of around 0.5 eV neutrinos become nonrelativistic very close to 
the epoch of matter-radiation equality. 

For all cases we provide in Table 1 limits on J2 after marginalizing over Aeflf and 
limits on Aeff after marginalizing over X^mj,. We stress that the neutrino mass scale 
can not be avoided as a standard cosmic fit parameter so that one should not derive 
limits on N^s while enforcing the neutrino masses to vanish. In practice, because 
and Aeflf are no longer degenerate, the allowed range for Aeff is not very different if one 
assumes = 0. 

The limits on A^^cfr differ little between our cases. Independently of the exact 
distribution of masses among the neutrino species we find 2.7 < N^e < 4.6 (95% CL), 
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except in Case 2 where the lower hmit is by definition 3 < A^eff • The bound is significantly 
stronger than the 2.5 < A^cfr < 5.6 found with WMAP-1 data and without inclusion of 
the BAO data [23]. The standard case of three massive neutrinos without modified 
number density and without additional radiation is well within the 95% CL range of 
Nqs, although it is just slightly outside the 68% CL range in Cases 1 and 2. Either way, 
the cosmological model with a nonstandard density of massive neutrinos or radiation is 
not significantly favored over the standard case. 
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